An Extended Abstract Submitted to 43 rd AIAA Thermophysics Conference, 

New Orleans, June 2012 


Electron Stark Broadening Database 
for Atomic N, O, and C Lines 


Yen Liu 1 2 , Winifred M. Hucf, Alan Wray 3 , and Duane F. Carbon 4 
NASA Ames Research Center, Moffett Field, CA, 94035, USA 


A database for efficiently computing the electron Stark broadening 
line widths for atomic N, O, and C lines is constructed. The line width 
is expressed in terms of the electron number density and electron- 
atom scattering cross sections based on the Baranger impact theory. 
The state-to-state cross sections are computed using the semiclassical 
approximation, in which the atom is treated quantum mechanically 
whereas the motion of the free electron follows a classical trajectory. 
These state-to-state cross sections are calculated based on newly 
compiled line lists. Each atomic line list consists of a careful merger of 
NIST, Vanderbilt, and TOPbase line datasets from wavelength 50 nm 
to 50 p,m covering the VUV to IR spectral regions. There are over 
10,000 lines in each atomic line list. The widths for each line are 
computed at 13 electron temperatures between 1,000 K - 50,000 K. A 
linear least squares method using a four-term fractional power series 
is then employed to obtain an analytical fit for each line-width 
variation as a function of the electron temperature. The maximum L 2 
error of the analytic fits for all lines in our line lists is about 5%. 


I. Introduction 

'T'HE emission and absorption coefficients in the radiative transfer equation for a 
A bound-bound transition depend on four quantities: the line position, the line oscillator 
strength or the Einstein coefficients, the line shape function, and energy state populations. 
The first two quantities are intrinsic properties of atoms and molecules. The NIST atomic 
spectra database [1], the Vanderbilt University (now also NIST) MCHF/MCDHF 
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database [2], and the Opacity Project [3] database TOPbase [4-5] provide the most up-to- 
date tabulation of these two quantities for many atomic species. However, their data are 
not all consistent to each other. For molecular species, on the other hand, most 
calculations still rely on approximate methods, such as the Dunham expansion for line 
positions and the R-centroid approximation for Einstein coefficients. The line shape 
function is determined by natural, Doppler, and collisional (pressure) broadening effects. 
While the natural and Doppler broadening shape functions are easy to determine, 
collisional broadening widths and shifts and rate coefficients to compute energy state 
populations require detailed quantum mechanical calculations of state-to-state cross 
sections involving the complete spectrum of energy states. Today, the most 
comprehensive source for collisional plasma (Stark) broadening is Griem’s book [6] 
published 37 years ago. His calculations were actually based on a limited number of 
states and covered only a small portion of transitions in the NIST line lists. In an effort to 
support the NASA Hypersonics Project’s technology challenges: (1) Develop new tools 
and technologies to enable high mass planetary entry and (2) Decrease uncertainty in 
aeroheating prediction by 50%, a new radiation code HyperRad is currently being 
developed at NASA Ames Research Center. Some of its key features are as follows: 

• Physics-based modeling for energy populations using coupled thennal, chemical, 
and radiative nonequilibria 

• Efficient and accurate linelist-driven database 

• Careful and complete merger of NIST, Vanderbilt, and TOPbase line datasets; the 
current database includes N, N + , O, 0 + , C, C + [7], and is easily extended to other 
atomic species 

• Molecular database based on ab initio calculations for electric dipole and 
quadrupole, magnetic dipole, and spin-forbidden transitions which include fine 
structure, pre-dissociation, and non-adiabatic corrections; the current database 
includes N 2 , 0 2 , C 2 , NO, CO, CN, N 2 + , CO + , C 3 , and C 2 H [8] 

• Microscopic collisional state-to-state cross sections currently based on a semi- 
classical method [9]; more accurate, quantal calculations are planned 

• Line broadening widths and shifts due to the Stark effect based on computed cross 
sections; Line -widths parameterized for efficiency up to 50,000K; maximum 
averaged error of parameterization is 5% 

• Nonequilibrium energy state populations based on a multi-group maximum 
entropy method; Macroscopic group equations and rate coefficients derived from 
moments of the master equations [10] 

• Three-dimensional, unstructured-grid approach for radiation transport applicable 
to any body shape and any CFD code 

• Designed ground-up for high-end computing in tenns of efficiency and 
parallelization 

In this paper, we present some HyperRad results for the computations of electron Stark 
broadening widths between 1,000K and 50,000K for atomic N, O, and C lines based on 
the Baranger impact theory [11-13] and a semi-classical approximation [14] using the 
complete set of energy levels and line lists compiled in Ref. 7. We also show that the 
Stark widths can be parameterized for efficiency using an analytical fit. 
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II. Impact Theory and Semi-classical Calculation 


The theory of collisional (pressure) broadening has been extensively studied by 
Baranger [11-13]. Using the impact approximation, an assumption that collisions take 
place sequentially and that two strong collisions never occur simultaneously, he showed 
that the collisional broadening shape function for a light-emitting particle (atom, 
molecule, or ion) in a gas of perturbers is a Lorentzian. The line shape of an isolated 
transition from an initial state i to a final state /can be expressed in terms of the half- 
width at half maximum (HWHM) y jf and frequency shift 6 jf as 

«,(»)-- / ; • (i) 

:t((o-(o if -6 if ) +y if 

Here 

(o v ~(o,-(o f (2) 

is positive for emission and negative for absorption. For a complete spectrum, one must 
sum over all possible final states and average over all initial states and all possible types 
of motion of the perturbers. The width y if can then be expressed in tenns of the perturber 

number density n p , the perturber velocity distribution function /(v) , and scattering cross 
sections o as 


Yif = y/oV (v ) (2^ + + °elJ dV • 

H i*f 


(3) 


Here o tj ( o fj ) is the inelastic cross section from the initial state i (final state / ) to a state 
j , and o elas is a combination of elastic cross sections a u and o ff , and the cross term. 


For electron-atom collisions, one of the most successful closed coupling approaches for 
computing the collisional cross section is the R-matrix with pseudo states (RMPS) 
method [15, 16]. However, the RMPS method can only handle a limited number of states, 
not the complete spectrum of states required in a Stark line shape calculation. We 
therefore adopt a semi-classical approximation [14] for our initial database construction. 
In the semi-classical approximation, the atom is treated quantum mechanically whereas 
the motion of the free electron follows a classical trajectory. The semi-classical inelastic 
cross section is thus given by 
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In these expressions, b designates the classical impact parameter of the perturber and 
( J,m ) is the total angular momentum and its projection on the axis normal to the 
collision plane. The symbol T denote the transition matrix operator 

{j\ T \ k ) = d jk ~ (j | exp[- l - V dt ] | k) , (6) 


where V is the interaction potential between the radiator and perturber. Collisional 
broadening in general is referred to as Stark broadening if the perturber is a charged 
particle. Here we consider electron Stark broadening. Let r and r' be position vectors of 
the free electron and atomic electrons, respectively. For electron-atom scattering, V is 
simply the Coulomb potential 


V(r,r') 



(V) 


where the free electron position r follows the classical trajectory. The computation of the 
cross sections in Eqs. (4) and (5) and details of differences from the previous theoretical 
calculations of Griem [6] and Riviere [17] will be presented in a companion paper [9]. 


III. Comparisons with Experimental and Theoretical Data 

For our calculations, we adopt the HyperRad atomic line lists for N, O, and C 
developed by Carbon and Huo [7]. The lists consist of a careful merger of NIST, 
Vanderbilt, and TOPbase line datasets from wavelength 50 nm to 50 urn covering the 
VUV to IR spectral regions. TOPbase lines were converted from LS coupling to .7.7 
coupling to include fine structure so that they are consistent with the NIST and 
Vanderbilt datasets. The line lists contain over 10,000 lines for each species in this 
spectral range. In Figure 1, we plot the Einstein A if coefficients for all the N lines below 5 

qm in the line list. In the figure, we also include the NEQAIR [18] line list for 
comparison. Note that NEQAIR has no data above 4 pm. 

In our calculations, we first compute the state-to-state cross sections in terms of 
electron velocity or energy, and the data are then convolved with the free electron 
distribution function to form the state-to-state rate coefficients. Using Eq. (3), we obtain 
the Stark width for a spectral line by first summing the rate coefficients for all other 
possible transitions from both the initial and final states of the line. The state-to-state 
cross sections and rate coefficients are also saved to compute the energy level 
populations. In order to evaluate the impact theory and the semi-classical approximation, 
we compare our calculations with several experimental data [19-21], as well as with 
Griem’s theoretical data [6], as shown in Figure 2 for several N lines in the UV region. 
Our calculations involve a complete set of energy states and therefore show relatively 
larger values than Griem’s data, but most are within the experimental uncertainties. 
Details of the comparisons, also including O and C lines, will be given in the final paper. 
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IV. Electron Stark Broadening Database 

As described in section II, the computation of the collisional broadening shape 
function involves very complicated and time-consuming procedures. Since the 
broadening mechanism is due to collisions, one must know the velocity distribution 
function of the perturber to compute the widths and shifts. In general, the distribution 
function can be expressed as a function of the macroscopic moments if moment equations 
of the Boltzmann equation are used to describe the fluid motion. Thus, for efficiency, the 
widths and shifts can be pre-computed and parameterized using macroscopic flow 
variables. In this paper, we employ the Maxwellian distribution characterized by an 
electron temperature T e to construct a database for Stark widths under the assumption 

that the free electrons are in translational equilibrium. The database for Stark shifts will 
be constructed in the future work. We first compute the widths at 13 T e positions, 

1,000K, 2,500K, 5,000K, 7,500K, 10,000K, 15,000K, 20,000K, 25,000K, 30,000K, 
35,000K, 40,000K, 45,000K, and 50,000K, for every transition line in our atomic N, O, 
and C line lists. Figure 3 shows the temperature variation of the Stark widths for all N 
lines. Each line width shown in the figure is normalized by its width at 10,000K. Due to 
lack of data, many computer codes employ a universal single temperature extrapolation 
formula for all lines to obtain widths at temperatures other than the reference one. For the 
purpose of comparison, such an extrapolation formula used in NEQAIR is also shown in 
the figure. Figure 3 also shows that the temperature variation of the Stark width is very 
smooth. This indicates that the line-width variation with T e can be parameterized for 

efficiency using either a spline or an analytical fit. The former is more accurate, but 
requires more storage for each line. In this paper, we construct an analytical fit for each 
line using the linear least squares method. The analytical fit is expressed in terms of basis 
functions (j>’"(T e ) as 

— K)-To?*"(t;) (8) 

n e m 

1 

For a four-term fractional power series (/)'" = (T )"' , m = 2 - 5 , the maximum L? error for 

all N lines is about 5%. Almost all large errors occur at temperatures below 2,500K, 
where Stark broadening is not important. At temperatures above 10,000K, the errors in 
general are below 1%. Figure 4 shows that the analytical fits essentially reproduce the 
semi-classical calculations at 12,500K for all N lines. Results for O and C lines will be 
presented in the final paper. 
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Figure 1. Atomic N lines 
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Figure 2. Comparison of Stark width with experimental and theoretical data. 
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Figure 3. Temperature variation of Stark widths for atomic N lines 



Figure 4. Comparison of semi-classical calculation and analytical fit at 12,500K. 
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